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ABSTRACT 

This  paper  is  concerned  with  asymptotic  expansions  of  solutions  of 
von  Kantian' s  swirling  flow  problem.  These  expansions  are  used  to  prove  the 
convergence  of  a  class  of  approximative  problems,  which  are  set  up  by 
substituting  for  the  infinite  interval  on  which  von  Kantian' s  problem  is  posed  by 
a  finite  but  large  one  and  by  imposing  supplementary  boundary  conditions  at 
the  far  end.  The  asymptotic  expansions  are  crucial  for  the  determination  of 
the  order  of  convergence.  Exponential  convergence  is  shown  for  well-posed 
approximative  problems. 

The  given  approach  is  applicable  to  general  autonomous  nonlinear  boundary 
value  problems  on  infinite  intervals,  for  which  the  von  Karman  problem  may  be 
considered  as  a  model  problem. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  is  concerned  with  the  von  Karman  swirling  flow  problem  which 
describes  the  velocity  field  of  a  fluid  over  an  infinite  rotating  disk.  It  is 
assumed  that  the  whole  half  space  over  the  disk  is  filled  with  fluid.  By 
restricting  the  class  of  admissible  solutions  von  Karman  (1921)  waf  able  to 
reduce  the  Navier-Stokes  equations  which  describe  this  fluid  dynamical 
configuration,  and  which  are  partial  differential  equations,  to  a  system  of 
ordinary  differential  equations  which  are  posed  on  an  infinite  interval.  This 
system  of  ordinary  differential  equations  can  well  be  considered  as  a  model 
problem  for  boundary  value  problems  on  infinite  intervals.  The  infinite 
interval  is  difficult  for  computation.  It  is  a  natural  procedure  to  restrict 
the  infinite  interval  to  a  finite  but  large  one  and  to  impose  additional 
conditions  at  the  far  end  which  should  describe  the  behavior  of  the  solution 
far  away  from  the  disk.  In  this  paper  information  about  the  asymptotic  (which 
means  far  away  from  the  disk)  behavior  is  obtained  and  used  for  the 
construction  of  these  far  away  boundary  conditions.  A  theorem,  which  implies 
convergence  of  the  'infinite'  solution  to  the  'infinite'  solution  as  the 
length  of  the  interval  becomes  infinite,  is  proved  and  exponential  convergence 
is  shown  for  suitable  approximating  problems.  The  extension  of  the  derived 
statements  to  more  general  problems  posed  on  infinite  intervals  is 
straightforward. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MBC,  and  not  with  the  author  of  this  report. 


ASYMPTOTIC  ANALYSIS  OF  von  KARMAN  FLOWS 
Peter  A.  Markovich 

1  .  Introduction 

von  Karman  showed  in  1921  that  the  Navier-Stokes  equations  for 
a  stationary  axisymmetric  flow  of  a  viscous  incompressible 
fluid  occupying  the  half  space  over  an  infinite  rotating  disk 
reduce  to  a  system  of  two  ordinary  differential  equations  when 
using  a  cylindrical  coordinate  system  (r,<p,z)  . 

The  disk  is  rotating  in  the  plane  z  =  0  around  the  z-axis  with 
angular  velocity  >  0  and  the  angular  velocity  of  the  fluid 
in  z  =  oo  is  yfiQ.  von  Karman  considered  only  the  case  y  =  0, 
but  the  generalization  to  y  £  JR  is  straightforward. 

With  the  normalization: 

1 

(1.1)  x  =  (nQ/v) 2  z 

where  v  is  the  specific  viscosity  of  the  fluid,  the  following 
similarity  equations  result  from  the  Navier-Stokes  equations: 

(1.2)  f’»  (x)  +  2f(x)f"(x)  =  (f’(x))2  -  g2(x)  +  y2] 

f  0  £  x  <  00 

(1.3)  g"(x)  +  2f(x)g'(x)  =  2f' (x)g(x)  J 
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The  velocity  field  of  the  fluid  in  the  (r,ip,z)  -  coordinate 
system  is  represented  by  (fiQrf  •  (x)  ,f2Qrg  (x)  ,-2  (vfiQ) 1  /2f  (x) )  . 
Appropriate  boundary  conditions  at  the  disk  are: 

(1.4)  f  (0)  =  a  £  3R,  f  *  (0)  =  0,  g(0)  =  1* 

A  non  zero  value  a  represents  blowing  from  the  disk  (a  <  O) 

or  suction  to  the  disk  (a  >  0) .  The  third  boundary  condition 
affecting  the  angular  velocity  implies  that  the  fluid  adhers 
to  the  disk.  Boundary  conditions  at  infinity  are 

(1  .5)  f'  («)  =  0,  g(=o)  =  y. 

For  y  =  0  McLeod  (1969a)  has  proven  the  existence  of  a  solu¬ 
tion  of  problem  (1.2),  (1.3),  (1.4),  (1.5)  for  all  at®. 

He  obtained  qualitative  properties  of  this  solution  (f,g) 
which  hold  independently  of  the  value  of  a: 

(1.6. a)  0  $  f'  <1,  (b)  lim  f(x)  =  f  €  3R 

OO 

X-*co 

(1.7)  g  >  0 ,  g '  <0 

The  solution  (f,g)  is  non-oscillatory ,  i.e.  neither  f'  nor 
g'  oscillate  around  zero  for  x  -»  <».  Moreover  f^  €  [-<®,oo) 
holds  for  every  solution  f  (see  McLeod  1970) ,  and  if  f^  6  IR+ 
then  f  -  f  ,f ' ,f " ,g,g'  decay  exponentially. 


For  y  *  0  McLeod  (1969a)  concluded  that 


(1.8)  signum  f'  (x)  \  const,  signum  g'  (x)  \  const  for  x  -»  » 

for  every  solution  (f ,g)  .  Furthermore  he  proved  the  existence 
of  a  solution  for  y  >  0  (1971)  independently  of  the  value  of 
a»  uniguenss  for  y  =  1  and  a  ^  0  and  a  nonexistence  theorem  for 
Y  =  -1  and  a  $  o.  Rogers  and  Lance  (1960)  heuristically  derived 
asymptotic  expansions  for  f  and  g  in  the  case  y  *  0  and  McLeod 

(1969b)  constructed  a  rig  -ous  proof  for  these  expansions 
which  imply  an  exponentially  dampened  oscillatory  decay  of 
f-f  and  g-Y-  But  his  proof  breaks  down  in  the  case  y  -  0. 

Lentini  and  Keller  (1980)  performed  a  numerical  study  of  the 
swirling  flow  problem  with  a  =  0  by  reducing  the  infinite  in¬ 
terval  [0,»]  to  a  finite  one  [0,X]  with  X»0  and  by  imposing 
two  additional  asymptotic  conditions,  at  x  =  X,  which  they  de¬ 
rived  from  the  requirement  that  unbounded  solution  components 

T 

of  the  in  ( f^O ,0 ,y ,0)  linearized  problem  vanish  at  infinity. 
Their  results  clearly  indicate  that  infinitely  many  solutions 
(f±,gi)  exist  for  Y  =  0  which  fulfill 

(1.9)  0  <  fi(<»)  <  f±_1  (~)  ,  i  =  2,3,... 

There  are  two  goals  of  this  paper.  The  first  is  to  derive  an 
asymptotic  expansion  for  f  and  g  in  the  case  y  =  0,  f  >  0, 
and  to  calculate  the  leading  coefficients  of  the  expansion 
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in  the  case  y  *  0.  The  second  goal  is  to  use  these  expansions 
in  order  to  prove  convergence  statements  for  reduced  problems 
posed  on  finite  intervals. 

The  organisation  of  the  paper  is  the  following:  in  paragraph 
2  we  reformulate  the  swirling  flow  problem  to  a  first  order 
boundary  value  problem  on  [0,«]  which  have  been  treated  by  de  Hoog 
and  Weiss  (1980)  andMarkowich  (1980)  and  collect  information  on 
the  solutions  of  these  problems,  in  paragraph  3  we  derive  the 
asymptotic  expansions  for  y  =  0,  in  paragraph  4  we  investigate 
the  case  y  *  0  and  in  paragraph  5  asymptotic  boundary  condi¬ 
tions  are  treated. 

The  method  employed  to  obtain  the  asymptotic  expansions  is 
linearization  of  (1.2),  (1.3)  around  (f  ,0,0,y,0)  (as  used 

OO 

by  McLeod  (1969))  and  a  contraction  mapping  theorem  which 
uses  the  closeness  of  the  solutions  of  the  linearized  and  non¬ 
linear  problems. 


2.  Reformulation  of  the  Swirling  Flow  Problem 


We  reduce 

(2.1) 


(1.2),  (1.3)  to  a  first  order  system  by  substituting 

y,  -  £,  y2  -  f.  y3  -  £",  y4  -  g,  y5  =  g’ 

y  *  (y, .y2.y3Iy4/y5)T 


and  get 


(2.2)  y'  =  F  (y ,  y) 


y3 

-2*1*3  +  y|  -  y*  +  Y2 
*5 

-2y1y5  +  2y2y4 


The  boundary  conditions  transform  to 


(2.3) 


'1 

0 

0 

0 

o' 

V 

0 

1 

0 

0 

0 

y  (0)  - 

0 

.0  0 

0 

1 

0. 

.1. 

(By  (0)  =  b) 


(2.4) 


lim 

X~*oo 


0  1  0  0  o' 

to] 

y  (x)  = 

o 

o 

o 

o 

y. 

(lim  Cy(x) 

X-KJO 


C(y)  )  . 


Because  of  (1.6.b)  and  (1.5): 

(2.5)  lim  y(x)  =  y  £  ro5  and  y  €  C([0,<»]) 

x-*« 


holds  and  (2.4)  can  be  changed  to  CY«=  c(y).  From  (2.5) 
and  (2.2)  we  deduce 


-5' 


r  1 


(2.6)  lim  y' (x)  =  F(y  ,y)  =  0 

OO 

X-*oo 

This  is  a  nonlinear  equation  which  can  be  solved  immediatly: 

(2.7)  y^  =  (fo,0,0,Y.O)T  with  €  3R. 

A  one  parameter  solution  manifold  y  =  y  (f  ,y)  which  is  con- 
sistent  with  (2.4)  is  specified  for  fixed  y  €  IR. 

The  next  step  is  to  calculate  the  matrix 

O  1  O  0  O 

OOIOO 
O  O  -2  f  -2y  O 

OO 

0  0  0  0  1 

0  2Y  O  0  -2  f 

OO 

which  has  the  eigenvalues  v^  defined  by  the  zeros  of  the  equa¬ 
tion 


(2.8)  A ( f  , y )  := 


3F(y«,(foo'Y)  »*> 

ay 


(2.9)  v[v2(v  +  2f  )2  +  4y2  ]  =  0 

OO 

These  eigenvalues  are; 


>- 

8 

4j 

a  0 

V2(fco'T) 

=  -  f 

OO 

+ 

C1(fco^) 

+  iC2(fco^) 

>- 

8 

44 

n 

> 

=  -  f 

OO 

+ 

ci  (fco'Y) 

-  ic^f^Y) 

’P 

u> 

(N 

1  > 

II 

'y- 

8 

44 

> 

=  -  f 

OO 

- 

o 

— fc 

8 

+  ic2(foo,Y) 

=  -  f 

OO 

- 

C^fco^) 

-  iVf^Y) 

=  Vfco'^ 
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I 


where  ,c2  are  defined  by 


(2.11) 


VVY)  -(i  (C*«T‘),/2*%/2 

«2‘^y)  =(l  (*i  +  *r')',I-$)U2 


If  y  =  0  all  eigenvalues  are  real  and  fulfill 


v2  >  0,  v3  >  0,  v,  =  v4  =  v5  =  0  if  <  0 

(2.12)  (a)  v1  h  v2  e  v3  =  0,  v4  <  0,  v5  <  0  if  f „  >  0 

V1  E  v2  E  V3  E  V4  E  v5  E  0  if  f~  =  0 


If  y  *  O  the  eigenvalues  v2'v3'v4'v5  have  a  nonvanishing  ima 
ginery  part  for  all  €  IR  and 


(2.12)  (b)  Re  v2  >  0,  Re  v3  >  0;  Re  v4  <0,  Re  v$  <  0  for  f^  e  3R. 
The  decrease  of  the  rank  of  A(f  ,y)  in  y  =  0  causes  the  dif- 

oo 

ferent  behavior  of  the  velocity  fields  for  y  =  0  and  y  *  0 
and  it  also  requires  different  treatment.  This  will  be  pointed 
out  in  the  paragraphs  3  and  4.  Now  we  transform  A(foo,T)  to  its 
Jordan  canonical  form  J(f  ,y) : 

(2.13)  Jtf^y)  =  E(fa5,Y)'1A(foo,Y)E(foo,Y) 
and  substitute 

(2.13)  Etf^Yjz  =  y  -  y^f^Y) 
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in  (2.2)  .  So  we  get  the  new  problem 


(2.14)  z’  =  J(f  ,y)z  +  h ( z , f  , y ) 


(2.15)  2  (<»)  =0 


where  h(2,f  ,y)  is  defined  by 


(2.16)  h(2,f  ,y)  =  E ( f  / Y )  F(E(f  ,y)z  +  y  (f  ,Y) /T)  -  J(f  ,Y)z 


Regarding  h(2,foo,Y)  for  fixed  y  and  f^  as  a  perturbation  to 
equation  (2.14)  we  come  up  with  an  ' inhomogenous 1  boundary 
value  problem. 


Using  a  suitable  norm  in  IR°  we  derive 


(2.17)  (a)  ||h(2,f  ,y)II  *  b.  (f  ,y)  II  z  ||  2 

OO  l  oo 


(2.17)  (b)  1 1  h  ( 2 1  ,  f  ^ ,  y  )  -  h(22 /f^Y)!!  S  b2  (f^Y)  (II  z^l  +  ll  z2ll  )  II  z1  -  z  2II 


where  b1  (f^Y)  and  b2(fco,Y)  are  independent  of  z  resp.  z 1  and 


In  order  to  solve  problem  (2.14)  we  have  to  find  an  appropriate 
particular  solution  and  therefore  let  us  briefly  discuss  the 


(general)  problem 


u'  =  Au  +  f(x),  x  ?  6  £  0 


(2.18) 


1 


where  the  Jordan  matrix  A  =  diag (Aq ,A+ ,A_)  is  partitioned  cor¬ 
respondingly  to  its  eigenvalues  will  real  part  equal,  larger 
and  smaller  zero. 

Pn  appropriate  choice  for  a  particular  solution  up  = 

T  T 

=(upo,up+,up_)  ,  where  f  =  (fQ,f+,f_)  is  partitioned  according 
to  A,  is  the  following: 

(a)  up  (x)  =  eAox  f  e  Aos  f  (s)ds 
o  J  o 

oo 

x 

(2.19)  (b)  up+(x)  =  eA+X  /  e“A+s  f+(s)ds 

OO 

(c)  up_(x)  =  eA-X  /  e  A-s  f_(s)ds 

6 

if  fQ  fulfills  f f  f o ( x )  ((  =  0(x  r  £)  where  c  >  O  and  r  is  the 
dimension  of  the  largest  Jordan  block  in  Aq.  This  choice  was 
suggested  by  Lentini  (1978)  and  de  Hoog  and  Weiss  (1980). 

In  operator  form  we  can  write 

(2.20)  up ( x )  =  (Hf) (x)  =  (HQfo,H+f+,H_f_) T(x) . 

H  is  a  linear  operator  on  the  space  of  all  function  f,  for 
which  f  €  C([8,°°])  and  II  f Q ( x ) II  =  0(x  r  £)  holds. 

If  A+,  A_  are  diagonalmatrices  the  following  estimates  are 
easily  derived. 
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(2.21) 


e  >  0 


(a)  II  (H+f+)  (x)ll  $  const  °||f+ll  ^ 

(b)  ||  (H  f  )  (x)ll  $  const  •  x_emax  ||tr+ef  (t)  II  , 

0  °  t^x  ° 

(c)  II  (H_f_)  (x) II  S  const  •  II  f _  II 

(d)  II  (H_f_)  (x)ll  S  const  -  || F_  1 1  x^x  v 

if  f_(x)  =  F_(x)x  v  with  y  >  0 

where  II  •  II  ^ a  denotes  the  max-norm  on  the  interval  [a,b]  . 

OO 

The  constants  depend  on  A  but  not  on  6  and  f.  If  f  e-A-sf(s)ds 

6 

exists  then  (H_f_) (x)  can  be  split  up  into: 

OO 

(2.22)  (H  f  ) (x)  =  eA-x  J  e“A-x  f  (s)ds  +  (H  f  ) (x) 

6 

where  H_  is  defined  by 

x  , 

(2.23)  (H_f  J (x)  =  eA-x  f  e  -  f_(s)ds. 

OO 

(Hf)  =  (HQfo ,H+f+,H_f_)  is  also  a  particular  solution  of 

(2.18).  if  f(x)  =  F(x)e  nox  with  1 1 F  ||  r ,  ,  $  const  and 

A_  =  -  n-j  •  I  with  n0/n.|  >0,  then  the  following  estimates  hold: 

(2.24)  II  (H  f  )  ( x ) II  $  const  II  F  II  r  ,e"nox,  n  >n, 

~  —  Lx  J  o  1 

(2.25)  ||  (H_f_)  (x)ll  ^  const  II  Fll  [6>xj  (e"n°x-e(ril“no)6e_niX) 

(2.26)  II  (Hf  )  (x)ll  $  const  II  F  II  r  ,e_nox 

o  o  Lx , °° J 
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After  straightforward  calculation  we  get 


(3.1) 


J(f  ,0)  = 


1 

0 


-2f 


-2  f 


(3.2) 


E(f  ,0)  = 


0 

1 

0 

0 

0 


0  0  1 

0  0  -2f 

0  0 

1  1  0 


4f2 

oo 


0  -2f  0 

OO 


(3.3) 


h(z,f  ,0}  = 


hi(z,foo,°) 


_  00 

-Sf2  Z5  (Zl 

+  z 

5)  +  (z2 

-2f~z5)2 

-  (z3 

+  z4)2)' 

_J_( 

2f  1 

.oo 

-8f»z5  (Z1 

+  z 

5)  +  (z2 

-2f»z5)2 

-  <z3 

+  z4)2) 

2T(4fcoZ4<Z1 

00 

fz5 

)  +2(z2 

“2f-*5) 

(z3  + 

Z4)  ) 

-2z4 

(Z1  +  z5)  - 

1 

f 

(z2-2f 

„  z5>  (z3  +z4) 

'2  Z5 

(z1  +z5)  + 

T 

4f2 

(  (z2  - 

2fcoZ5)2' 

(z3  + 

Z4)  2) 

) 


T 

assuming  f^  >  0.  Here  z  =  (z1  ,z2,z3,z4,z5)  holds.  The  general 
solution  of  problem  (2.14)  for  fixed  f^  is: 


(3.4) 


z  =  exp  ( J  (f^OJx) 


+  (Hh(z,f  00,0) )  (x) 


The  proposition  (2.15)  holds  iff  vi  =  v2  =  v3  =  0  because  2 
decays  exponentially  and  the  estimates  (2.17) (a)  and  (2.21)  (b) 
hold. 


For  f^  €  3R  ;C1»C2  e  3R  fixed  we  define 


(3.5) 


(*o(z'fJ  >  (x)  “ 


0 

0 


0 

0 


0  0 
—2  f<»x 


0  e 


0 

-2f0Ox 


(II) 


+  (Hh ( z , f  ,0) ) (x) 


as  an  operator  on  the  Banachspace 


(3.6)  (A-  =  { u  |  u  ( x)  =  U(x)e"2fooX,  U€  C.  ([6,oo))  },  Hull  =  Hull  ) 

foo  b  [6,co] 


where  C^US,*®))  is  the  space  of  all  functions  which  are  conti¬ 
nuous  on  [  6,oo)  and  bounded  on  [  6 ,<*» ]  .  6  ?  0  will  be  chosen  sui¬ 
tably. 


ip  (  • , f  )  maps  Ac  into  Ar  because  of  the  estimates  (2.17)  (a), 

O  oo  *  z  L 

oo  oo 

(2.21) (a),  (b)  and  (2.25).  Every  solution  z  is  a  fixed  point 
of  ,J'o(,»fO0)  and  the  existence  and  uniqueness  of  this  fixed  point 
is  shown  by  the  contraction  mapping  theorem,  which  is  appli¬ 
cable  for  6  sufficiently  large  on  any  sphere  with  center 
(0,0,0,e_2f°°x51  ,e"2fooXC2)T  because  the  estimates  (2.17)  (b), 
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(2 .21 )  (a)  , (b)  and  (2.25)  hold.  As  mentioned  before  all  so¬ 
lutions  z  decay  exponentially,  however  we  have  to  prove  that 
there  is  no  fixed  point  of  (  * ,  f^)  which  decaysslower  than 
e~  2  fcoX . 

This  is  done  by  regarding  <|i  (*,f  )  as  an  operator  on  the  Ba- 
nach-space 

(3.7)  (B£  =  (vlv(x)  =  V(x)x'2_e,V  e  Cb([  6,~)  )  ,e  >  0},  II  v  II  *=IIVll[6 

The  contraction  mapping  theorem  can  be  applied  again  for  6 
sufficiently  large  and  because  <=  B£  holds  for  every 

oo 

t  >  0 , f  >  0  the  established  solutions  z  6  Ax  are  unique  in  B„  . 

oo  r  t 

oo 

The  only  difference  in  the  proof  of  the  contraction  properties 
of  (*,f  )  in  B  and  A,.  is  that  the  estimate  (2.21)  (d)  has  to 

O  8  T. 

oo 

be  used  instead  of  (2.25). 

Assume  now  that  f  <  0.  If  z  E  B„  then  Hh(z,f  ,0)  £  Bo  and 

oo  oo  £ 

therefore  v1=v2=v2=^i=^l=0  has  to  be  and 

by  the  contraction  mapping  theorem  z  =  0  is  the  only  fixed  point. 
For  =  0  the  same  argument  is  valid  in  B-|+ef  because  the 
largest  Jordanblock  of  J(0,0)  is  three  dimensional. 

We  exclude  f  <0  because  we  only  look  for  solutions  in  B„  for 
f  *  0  and  B. ,  for  f  =0. 

oo  1+8  oo 

—  2  f  x 

In  order  to  investigate  whether z might  decay  faster  than  e  00 

we  substitute  H_ ,  which  is  defined  by  (2.23),  by  H_  and 

H  =  (HQ,H+,ff_)  for  H.  Because  of  the  estimate  (2.24)  with  nQ  =  4^ 

and  n.  =  2f  and  because  of  (2.2b)  we  conclude  that 
1  oo 


(3.8)  ||{Hh(z,f  ,0) )  (x)||  =  0 (e”4f°°x)  for  z  €  A- 

oo  z 

i 

holds.  We  have  cut  off  the  0(e"’2^0°x)  term  in  H_. 

The  (modified)  operator  ( * » f^)  is  defined  by 

0  0 

0  0 

(3.9)  (*_(z,f  ))(x)  =00 

e“2fooX  Q 

0  e~2 

The  contraction  argument  in  Af  applies  again  and  because  a 

CO 

fixed  point  z  of  (  • , f^)  establishes  a  solution  of  (2.14), 
(2.15)  the  equivalence  of  the  solution  manifolds  defined  by 
the  fixed  points  of  and  ^  for  £  3R+  and  U.^,^)  £  IR2. 
follows  because  z  €  Af  and  because  (2.22)  holds. 

•*-  OO 

Assume  now  that  (w1 ,w2)  =  (0,0).  Then  the  unique  fixed  point 
of  ‘^O(*,foo)  is  z  =  0  implying  the  trivial  solution  f  =  f^, 
g  s  0. 

Moreover  assume  that  w^  =  0. 

Then  start  the  iteration  zn+^  =  $  ( zn,f  )  with  some  vector 

O  oo 

z°  =  (z°,  z^  ,0,0,  Zj) T  £  Af  ,  (3.3)  implies  that  h(z°,foo,0)  = 

OO 

=  (z1 ,z2,0,0,z5)  . 

Therefore  the  (unique)  fixed  point  z*  of  ^0(*,f00)  fulfills 

z*  =  lim  zn  =  (z* ,  zSj  ,0,0,  zjj!) T  £  Af  ,  which  implies  after  pre- 
n-w°  “ 

multiplication  with  £(^,0)  that  g  =  0. 
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From  the  definitions  of  and  h  we  derive 


(a)  Zl  =  C(e"4foaX) 

(b)  z 2  =  0(e'4f“x) 

(3.10)  (c)  z 3  =  0(e"6f»x) 

(d)  z4  =  e~2f°°x  w  1  +  0(e”6f°°x) 

<•>  *5  *  e'2f“*  w2  +  (jr  wi  +  5pr  wj)  e'4*"*  +  ofe-6^”*) 

v  oo  oo  7 

_  4  f  V 

(3.10)  (c)  and  (d)  hold  because  the  terms  of  order  e  can- 

cel  in  h, (z,f  ,0)  and  h.(z,f  ,0). 

J  OO  4  oo 

By  applying  the  identity  y  =  Eff^OJz  +  yoo(fot)fO)  •  integrating 
2 

f"'(x)  =  4^  z,-(x)  twice  and  by  continuation  from  [  6,°°]  to 
[0,oo]  we  get 


THEOREM  I 

Every  solution  of  (1.2),  (1.3),  (1.5)  with  f^  >  0 

and  y  =  0  fulfills: 


f(x)  =  f  +  e 


-2f~x 


W2  +  {sf  W2  + 32 f7  Wl) 

'  OO  OO  • 


-4f  X 
e  oo  + 


+  0(e"6fooX)  ,  x 


g (x)  =  e_2f°°x  w1  +  0(e'6fa>x),  x 


where  w^ ,  €  IR  .  w^  =0  implies  g  e  o,  w^  =  w^  =  0 

implies  f  =  f  ,  g  =  0.  Therefore  no  nontrivial  so¬ 
lution  (f,g)  oscillates  for  large  x,  the  convergence 
of  f  -  f  and  g  to  0  for  x  -*  00  is  monotonic. 

OO  3 
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The  boundary  value  problem  {2.2),  (2.3),  (2.4) 
regarding  the  number  of  boundary  conditions  at 
x  =  oo,  because  the  three  parameters  (foo,w1,w2) 
termined  from  the  three  boundary  conditions  at 
(2.3)  . 


is  well-posed 
x  =  0  and 
have  to  be  de- 
x  =  0  given  by 
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In  this  case  the  eigenvalues  of  A(f  ,y)  are  distinct  and  we  get 


(4.1)  Jtf^Y)  =  diag(0,v,(f  .y) ,v,(f  ,T) ,v.(f  ,T) ,v,(f  ,y)) 


1 

0 

(4.2)  E ( f  ,y)  =  0 

oo 

0 

0 


1 

1 

1 

1 

V2 

V3 

V4 

V5 

2 

V2 

V2V5 

V? 

V3V4 

< 

v|v2 

2y 

V2V5 

2y 

V3V4 

2y 

V4V3 

2Y 

v|v2 

2y 

2Y 

2y 

2Y 

Pursuing  as  in  paragraph  3  and  using  that  v3  =  v4  we  get 

0  O 

O  0 

O  0 

ev4(fco,Y)x  Q 

0  ev4(foo,Y)x 

Defining  ♦  (z^f^)  by  the  right  hand  side  of  (4.3)  as  an  opera¬ 
tor  on  the  Banachspace : 

(4.4)  (Af  t  =  {u|u(x)«U(x)e(ReV4  (foo,r)  )x,U€Cb([  6,®)  } ,  ||  u  II  =  II U  II  [  6f0(>]) 
we  conclude  as  we  did  in  paragraph  3  that  ip  (-,f  )  is  a  contrac- 

T  _°° 

tion  on  any  sphere  with  center  (0,0,0,eVl*x51 ,  eVl*x5i;)T,  if  6  is 
sufficiently  large. 
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(4.3)  z  « 


,*? 


1 

I 

i 

McLeod  (1969b)  proved  that  there  is  no  solution  (f,g)  that  decays 
slower  to  (f^Y)  than  eReV3X.  We  want  to  prove  now  that  there 
is  no  solution  z  of  (4.3),  which  decays  faster  than  that.  There- 

1 

fore  we  again  substitute  H  for  H  and  as  in  paragraph  3  we  find 

i 

l  that  w1  =1*2  =  0  yields  z  =  0  as  the  unique  fixed  point  of 

(4.4)  (?Y(z,fJ)(x)  = 

So  we  have  established  the  existence  of  a  unique  (f^w^  ,w2) - 

2 

solution  manifold  for  f  €  1R  ,  (w, ,w„)  £  I  . 

oo  I  Z 

The  real  solution  manifold  y  =  E(f  ,y)z  +  y  (f  ,y)  is  obtained 
by  choosing  w1  =  ^(0.^102)  and  w2  =  w1  .  The  three  parameters 
are  now  (f  ,0,,0_)  e  IR3  .  In  order  to  show  this,  we  set 
u  =  E(f  ,y)z  and  get  a  fixed-point  equation  for  u  by 

I 

i  (4.5)  u  =  E(f  ,y)?v(E_1  (f  ,Y)vi,f  ) 

OO  Y  OO  OO 


o  0 

0  0 

0  0 

5V4X  0 

n  3X 


/W  \ 

^w2j  +  (ffh(z,foo,Y) )  (x) 


and  therefore  u  =  lim  u  where 

n-*«> 

(4.6)  un+1  =exp(A(f  ,Y)x)E(f  ,y) 


r0  ' 

0 

0 

w. 


+  E(foo,Y)  (Hh  (E~ 1  (foo,Y)un,foo,T))  (X) 


holds.  The  first  term  on  the  right  hand  side  is  real  because 

w.  =  w„ ,  the  second  is  real  if  u  is  real  because 
I  e.  n 
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(4.7)  E  (foo,Y)exp(J(feo,r)  (x  -  s)  )E-1  (f^Y^xpfAff^Y)  (x  -  s) ) 

holds.  This  yields  inductively  that  u  as  well  as  y  =  u  +  yB 
are  real,  if  the  starting  function  u°  is  real.  Therefore  we 
get 

THEOREM  II 

Every  solution  of  (1.2),  (1.3),  (1.5)  with  y  /  0 

2 

fulfills  for  some  ( 3 ^ , 3 2 ^  £ 

f(x)  =  f w  +  e"(f«+ci(f“'Y))x(01  cos  (c2(foo,Y)x)  + 

+  02  sin  (c2(foo,Y)x))  +  0(e“2  (f“+Cl  (f“'Y)  )x)»  x  +  » 
g  (x)  =y+2^t  e'(foo+Cl  (foo'Y))x(ai  cos  (c2(foo,Y)x)  + 

+  a2  sin  (c2(foo,y)x))  +  0(e-2  (foo+c2  (foo'Y) )  x)f  x  -  » 

where  a,  =  (f^Y)  =  Re(v2  (fa>,Y)v3(foo,Y)  (3  1-i32) } 

and 

a2  =  a2{f«>'y)  “  “  Im(v4(foo'T)v3(foo,Y)  (B1  -  i02)  )  . 

Furthermore  01  =  32  =0  yields  f  h  f^  and  g  =  y . 

All  nontrivial  solutions  f  and  g  oscillate  for 
large  x  in  the  sense  that  there  are  sequences 
x .  ,  x .  €  IR  with  lim  x .  =  lim  x ,  =  + «.  which  are 

1  1  i-wo  1  i-*oo  i 

the  only  zeros  of  f  -  f^  resp.  g-Y  and 
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xi+1  "  xi 
[xi+1  "  xi ) 


Wf  +  i 


holds  for  i  ■*  ®. 


Again  the  boundary  value  problem  on  [0,°°]  is  well  posed  re¬ 
garding  the  number  of  boundary  conditions  at  x  =  0  and  x  =  «>. 

A  convergent  expansion  for  f  and  g  of  the  form 

OO 

I  c  (f  ,w.  ,w„ ) exp (-n  (f  +c.{f  ,y))x)  can  be  obtained  for  all 
n  '  OOr  12  OO  \  CO*  ' 

n-o  _  i  m  ^  —  q 

T  £  1  from  the  iteration  z  =  ^  (z  ,f  )  by  choosing  z  =0 

T  00 

as  starting  function. 
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5.  Asymptotic  Boundary  Conditions 

Boundary  value  problems  on  inifinite  intervals  are  often  sol¬ 
ved  by  restricting  to  a  finite  but  large  interval  and  by  im¬ 
posing  supplementary  boundary  conditions  on  the  right  end  (see 
Lentini  (1978),  de  Hoog  and  Weiss  (1980-)  and  Markowich  (1980)). 
The  approximative  problems  for  the  swirling  flow  problem  have 
the  following  form: 

(5.  1 )  u'  =  F(u,y) 

(5.2)  Bu (0)  =  b  0  $  x  $  X,  X>>0 

(5.3)  S(u(X),Y)  =  0,  S:3R5  -*  ]R2  for  fixed  Y  £  ® 

The  two  additional  asymptotic  boundary  conditions  (5. 3)  shall 
express  the  desired  behaviour  of  the  exact  solution  y  for 
large  x  so  that  they  shall  assure  convergence  of  u  to  y  in 
the  following  sense: 

(5.4)  Ilu(  -,X)-yll[0  -»  0  as  X  -*  °°. 

Moreover  the  order  of  convergence  should  be  reasonably  large 
in  order  to  enable  us  to  compute  the  solution  u(x,X)  of  the 
two  point  boundary  problem  (5.1),  (5.2),  (5.3)  sufficiently 
fast. 

From  Spijker's  nonlinear  stability-consistency  concept  (see 
Spijker  (1972),  paragraph  3.1)  we  conclude  that  it  is  suffi¬ 
cient  for  the  convergence  to  require  that  1)  and  2)  hold,  where: 
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1)  The  problem 


(5.  5) 


v  * 


3F 

3y 


(y  (x)  #y)  v  +  f  (x) 


0  ^  x  ^  X 


(5.  6)  Bv  =  0 

(5.  7)  f|  (y(X)  ,Y)V(X)  =  Y(X) 

is  for  all  f  e  C([0,x]),  Bem3,Y(X)  e  m2  uniquely  soluble  and 
the  solution  v  fulfills  a  stability  estimate  of  the  form: 

(5.8)  il  vll  [0,x]  *  0,3  11  +xr'1|,Y(X)  11  +  X2r  1||f"[0,X]) 

where  r  is  the  dimension  of  the  largest  Jordanblock  belonging 
to  zero  eigenvalues  of  Atf^Y)  ,  that  means  r  =  1  for  y  *  0,^  £  M 
and  r  =  2  for  y  =  0,  f  >0. 

’  CO 


2)  the  consistency  assumption 

(5.9)  II  S  (y  (X)  ,Y)II  =  0(X"2  (2r_1)_e)  ,  e  >  O 

holds. 

At  least  convergence  of  the  order 
(5.  10)  X2r_1l|S(y(X)  , y ) II 


results  and  the  approximative  problems  (5.1),  (5.2),  (5.3)  are 
for  all  sufficiently  large  X  soluble.  This  solution  u(x,X) 
is  unique  in  a  ball,  whose  center  is  the  restriction  of  y  to 


[0 ,  X]  and  whose  radius  is  smaller  than  const.  X  ^2r-1^  (in  the 

space  (  [0/ X] )  with  Hull  =  Hull  ^  +  llu'll  )  where  const, 

is  sufficiently  small. 

The  basic  assumption  is  to  require  the  solution  of  (2.2),  (2.3), 
(2.  4)  to  be  isolated,  i.  e.  the  linearized  problem 

(5.  11)  w'  =  —  (y(x) ,y)w,  0  ^  x  $  ®,  w  €  C([0,»]) 

°y 

(5.  12  )  Bw (0)  =  0 

(5.  12)  Cw(oo)  =  o 

has  the  unique  solution  w  =  0,  so  bifurcation  and  limit  points 
are  excluded  from  the  convergence  analysis.  Writing  (5.11)  as 

W'  =  3y  (yoo'rJw  +  G(x)w 
with  G(x)  =  0(||y(x)  -  y^Cf^Y)  II  we  get 

(5.  13)  W  =  E(foo,Y)exp(J(foo,Y)x)C  +  E  ( f ^ ,  Y )  (HE-1  (f  ,,y)Gw)  (x) 

where  £;  =  U  y  ^ 3f  $4  >  •  By  applying  (5.12)  we  conclude 

that  ?2  =  C3  =  0  for  y  =  0  as  well  as  y  *  0. 

The  general  solution  of  (5.11),  (5.12)  is 

'  1  0  0] 

0  0  0  n1 

(5.14)  w(x)  =  E  ( f  , y )  4>  ( x )  0  0  0  n- 

0  10^ 

.0  0  lJtn3 

where  $  (x)  is  a  fundamental  matrix  of 


U'  =  E-1  (f00/T)|^(y  (x)  ,Y)E(foo,Y)u 

Therefore  y  is  isolated,  iff  the  3x3  Matrix 

‘1 

0  0 

(5.15)  BE  ( f  ,y)$(0)  0 

oo 

l°0 

is  regular. 

From  this  proposition  and  certain  assumptions  on  (y(X) ,y) 
Markowich  (1980)  concluded  the  unique  solvability  of  (5.5), 
(5.6),  (5.7)  for  all  (f,&,Y(x))  6  C([0,xl)  x  and  the  esti¬ 
mate  (5.8)  similarly  to  de  Hoog  and  Weiss  (1980)  who  dealt  with 
problems,  where  F  (y  )  has  no  eigenvalue  with  real  part  zero. 
These  assumptions  are: 

0  O' 

(5.16)  llf||(y(x) /Y)E(f  ,Y)  0  1  )_1||  =0(1)  for  X  -  - 

\<*y  00  00' 

.0  0. 

(5.17)  II  ||(y(X)  ,Y)  II  =0(1)  for  X-*- 

'  1 ' 

8S  ®  — r 

(5.18)  11  (y ( X)  ,Y)E(f  ,y)  0  II  =  0 (X  r)  for  X  -  » 

ay  0 

.0. 

Because  y(X)  decays  exponentially  it  is  sufficient  to  require 
~(y  (f  ,y)  ,Y)  to  fulfill  (5.16),  (5.17),  (5.18)  instead  of 

dy  ®  00 

!§(y(X) ,Y) .  Obviously  (5.9)  is  fulfilled,  if 

dy 


\ 


holds,  and  then  the  order  of  convergence  is  at  least 
X»exp (- (f  +  c..  (f  , y)  ) X)  for  y  *  0  and  X"^exp(-2f  X)  for  y  =  0 
because  of  the  expansions  for  y.  So  (5.16),  (5.17),  (5.18)  with 
y^f^Y)  substituted  for  y(X)  and  (5.19)  specify  a  class  of  ad¬ 
missible  supplementary  boundary  conditions.  A  natural  choice  is 

S(U(X),T)  =  [g  o  0  °  o]U<X>  '  (?)  *  Cu(X>  '  =  <*> 

Obviously  (5.19)  is  fulfilled  and  it  is  checked  by  a  straight¬ 
forward  calculation  that  (5.16),  (5.17),  (5.18)  hold  for  Y  =  0 
as  well  as  y  *  0.  From  the  asymptotic  expansion  given  in  Theo¬ 
rems  1  and  2  we  conclude  that  the  order  of  convergence  is  mini¬ 
mal  for  this  choice,  i.e.  Xexp(-(f  +c.(f  ,y))X)  for  y  *  0  and 

3 

X  exp(-2fooX)  for  y  =  0,  but  even  these  linear  inhomogenous  supp¬ 
lementary  boundary  condition  yield  a  reasonable  -  exponential  - 
order  of  convergence . 

Lentini  and  Keller  (1980)  used  the  following  asymptotic  boun¬ 
dary  conditions 


Sp  (u  (X)  / T) 


(u1 (X)+cl (u1 (X) , y) )u2(X)+u3(X)-c2(u1 (X) , y) (u4(X)-y)1 
c2(u1 (X) ,y)u2(X)+(u1 (X)+c1 (u1 (X) ,y) ) (u4(X)-y)+u5(X), 


for  y  *  0  and  their  limiting  forms 


Sp (u (X)  ,0) 


2u1  (X)u2(X)  +u3(X)' 
,2u1  (X)u4(X)  +u5(X), 


0 
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for  Y  -  0.  These  so  called  projection  conditions  are  derived  by 
setting  the  zero  these  solution  components  of  the  linear  prob¬ 
lem  (2.14)  (without  the  1 inhomogenous '  term  h (z , f^, y) )  and  (2.15) 
which  do  not  automatically  vanish  at  infinity.  Checking  the  con¬ 
ditions  (5.16),  (5.17),  (5.18)  and  (5.19)  yields  the  stability 
and  consistency  of  the  projection  conditions. 

In  order  the  determine  the  exact  order  of  convergence  which  is 
given  by  (5.9)  we  substitute  the  asymptotic  expansions  for  y 
from  Theorem  1  in  the  case  y  =  0  and  from  Theorem  2  for  y  *  0 
for  u(X)  into  the  projection  conditions  and  get 

(5.20)  HSp  (y  (X)  ,Y)||  =  0(exp(-2  (fw  +  c.,  (f^y)  )X))  for  y  £  1 

Therefore  the  projection  conditions  cancel  the  first  two  terms 
of  the  series  expansion  for  y.  The  expression  (5.20)  also  ex¬ 
plains,  why  Lentini  and  Keller  (1980),  who  assumed  no  blowing 
or  suction  at  the  disk  had  to  increase  the  length  of  the 

interval  X  in  order  to  achieve  the  same  accuracy  when  computing 
higher  order  Karman  swirling  flows.  Because  of  (1.9)  the  se¬ 
quence  fi(»)  =  f*  for  i  *  1,2,...,  is  decreasing,  therefore  the 
order  of  convergence  decreases  too  unless  X  is  increased. 


Acknowledgement . 


I  am  very  grateful  to  Professor  Richard  Weiss  from  the  In 
stitute  for  Numerical  Mathematics  at  the  Technical  univer 
sity  of  Vienna  for  many  hours  of  stimulating  discussions. 


References 


F.R.  de  Hoog  and  R.  Weiss  (1980) 

An  Approximation  Method  for  Boundary 
Value  Problems  on  infinite  Intervals. 
Computing  24,  pp.  227-239. 

T.  von  Karman  (1921)  Uber  laminare  Reibung. 

Z.  Angew.  Math.  Mech.  1,  pp.  232-252. 

M.  Lentini  (1978)  Boundary  Value  Problems  over  Semi- 

finite  Intervals. 

Ph.D.  Thesis,  California  Institute  of 
Technology,  Pasadena,  CA  91125. 

M.  Lentini  and  H.B.  Keller  (1980) 

The  von  Karman  Swilling  Flows. 

SIAM  J.  Appl.  Math.,  Vol.  38,  No  1., 
pp.  52-64. 

P.A.  Markowich  (1980)  Randwertprobleme  auf  unendlichen  In- 

tervallen. 

Dissertation,  Technische  Universitat 
Wien. 

J.B.  McLeod  (1969a)  Von  Karman's  Swirling  Flow  Problem. 

Arch.  Rat.  Mech.  and  Anal.,  Vol  33 
pp.  91-102. 

J.B.  McLeod  (1969b)  The  asymptotic  Forms  of  Solutions  of 

von  Karman's  Swirling  Flow  Problem. 
Quart.  J.  Math.  Oxford,  (2)  20, 
pp.  483-496. 

J.B.  McLeod  (1970)  A  Note  on  rotationally  symmetric  Flow 

above  an  infinite  rotating  Disk. 
Mathematika,  Vol.  17,  pp.  243-249. 

J.B.  McLeod  (1971)  The  Existence  of  axially  symmetric 

Flow  above  a  rotating  Disk. 

Proc.  Roy.  Soc.  A  324,  pp.  391-414. 

M.H .  Rogers  and  G.N.  Lance  (1960) 

The  rotationally  symmetric  Flow  of  a 
viscous  Fluid  in  the  Presence  of  an 
infinite  rotating  Disk. 

I.  Fliud  Mech.  7,  pp.  617-631. 

M.N .  Spijker  (1972)  Equivalence  Theorems  for  nonlinear 

finite  Difference  Methods,  in: 
Numerische  Lfisung  nichtlinearer  par- 
tieller  Differential  und  Integrodif- 
f erentialgleichungen . 

Lecture  Notes  in  Mathematics,  Vol.  267. 


-28- 


SECURITY  CLASSIFICATION  Or  THIS  RAGE  r,1>x*i  J>*»a  Kntned) 

I  REPORT  DOCUMENTATION  PAGE 


t.  REPORT  NUMBER 

2156  / 


[«.  TITLE  (and  Subtitle) 


2.  GOVT  ACCESSION  NO 


ASYMPTOTIC  ANALYSIS  OF  von  K ARMAIS  FLOWS, 


7.  AUTHOR(i) 

)  A  Peter  A. /Markovich 
1  .  1 

9-  PERFORMING  ORGANIZATION  NAME  AND  ADp R E S S 

Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison,  Wisconsin  53706  _ 

It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


2=1 -Lk£kLlt _ 

K1£AD  INSTRUCTIONS 

_ Df-PORK  COUPI.r.TlNT,  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

Summary  Report  -  no  specific 
_ reporting  period _ 

6.  PERFORMING  ORG.  REPORT  NUMBER 
a— -CONTRACT  OR  GRANT  NUMBER^iJ 

/£/  * 

DAAG29-80-C-0041* 

£-MC^  9-27062 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

1  -  Applied  Analysis 


See  Item  18  below 


lYW  //! 


12.  REPORT  OATE 


)  Dec  mfg&x  1980 

13.  NUMBER  OF  PAGES 
28 


14.  MONITORING  '>GENCY  NAME  ft  ADDRESSfif  di  tie  rent  from  Controlling  0///ceJ  |  15-  SECURITY  CLASS,  (of  thie  report) 


^  l  t  •  I  £0  h  '  r 

_ _  ftu  _ 

16.  DISTRIBUTION  STATEMENT  (at  this  Report)  ' 

Approved  for  public  release;  distribution  unlimited. 


UNCLASSIFIED 


15a.  DECLASSIFICATION  OOUNGRADir  G 
SCHEDULE 


[T7.  DISTRIBUTION  STATEMENT  (at  fat  ab.lracf  entered  In  Block  30,  II  dlllerent  from  Report) 


18.  SUPPLEMENTARY  NOTES  ,  ,  . 

U.S.  A.  itrmy  research  Office 
P.O.  Box  12211 
Research  Triangle  Park 
Worth  Carolina  27709 


National  Science  Foundation 
Washington,  D.  C.  20550 


19.  KEY  WORDS  (Continue  or  reverie  wide  if  necenery  end  Identify  by  block  njoiber) 

Navier  Stokes  equation,  singular  points,  asymptotic  properties, 
asymptotic  approximation 


"JlT  ABSTRACT  (Continue  on  reverie  tide  It  necesimry  end  Identify  by  block  nu7ib«rj 

>^rhis  paper  is  concerned  with  asymptotic  expansions  of  solutions  of  von  Kantian1  s 
swirling  flow  problem.  These  expansions  are  used  to  prove  the  convergence  of  a 
class  of  approximative  problems,  which  are  set  up  by  substituting  the  infinite 
interval  on  which  von  Kantian’ s  problem  is  posed,  by  a  finite  but  large  one  and  by 
imposing  supplementary  boundary  conditions  at  the  far  end.  The  asymptotic  expan¬ 
sions  are  crucial  for  the  determination  of  the  order  of  convergence.  Exponential 
convergence  is  shown  for  well-posed  approximative  problems. - "(continued)  - 


DD  ,^M73  1  473  EDITION  OF  I  NOV  «»  IS  OBSOLETE 


UNCLASSI FIED  ^  J,  (J  Q  V, 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (HJian  Data  Entered) 


20.  (Abstract  continued) 

The  given  approach  is  applicable  to  general  autonomous  nonlinear  boundary 
value  problems  on  infinite  intervals,  for  which  the  von  Karman  problem 
may  be  considered  as  a  model  problem. 


